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We combined quantitative relaxation rate (/?, = 1/7", ) mapping — to 
measure local myelination — with fMRI-based retinotopy. Gray- 
white and pial surfaces were reconstructed and used to sample 
at different cortical depths. Like myelination, /?, decreased from 
deeper to superficial layers. /?, decreased passing from VI and MT, 
to immediately surrounding areas, then to the angular gyrus. High 
/?, was correlated across the cortex with convex local curvature so 
the data was first "de-curved". By overlaying R : and retinotopic 
maps, we found that many visual area borders were associated 
with significant /f, increases including VI, V3A, MT, V6, V6A, V8/ 
V01, FST, and VIR Surprisingly, retinotopic MT occupied only the 
posterior portion of an oval-shaped lateral occipital /?, maximum. 
/?, maps were reproducible within individuals and comparable 
between subjects without intensity normalization, enabling multi- 
center studies of development, aging, and disease progression, and 
structure/function mapping in other modalities. 

Keywords: MT, myelination, parcellation, surface reconstruction, visual 
areas 



Introduction 

Cortical areas are best defined and recognized on the basis of 
multiple converging techniques (Allman and Kaas 1971). Five 
measures developed during invasive experiments on animals 
include 1) receptotopic organization (e.g. retinotopy), 2) archi- 
tectonic features, 3) connection patterns, 4) neurophysiological 
properties, and 5) effects of localized lesions. These measures 
have each been adapted to human brains. Studies on the effects 
of brain damage in humans have a long pedigree, recently sup- 
plemented by transcranial magnetic stimulation studies in 
normal subjects. Cortical-surface-based functional magnetic res- 
onance imaging retinotopy (Engel et al. 1994; Sereno et al. 1995; 
DeYoe et al. 1996) is now established and often combined with 
functional studies, which constitute the overwhelming majority 
of in vivo human neuroimaging experiments. Connection pat- 
terns are beginning to be addressed in vivo with diffusion- 
imaging-based fiber-tracking methods. But despite the historical 
precedence of postmortem human architectonic studies, in vivo 
human architectonics has perhaps been the least well-adapted 
of the 5 measures. 

Postmortem cortical parcellation remains an difficult 
problem a century after the publication of Brodmann's classic 
cytoarchitectonic maps of the cortex in humans and other 
species. Related early work using cytoarchitectonic and mye- 
loarchitectonic methods soon resulted in a plethora of human 



cortical area maps (reviews: Zilles and Amunts 2010; Geyer 
et al. 2011). But human cortical parcellation then went into 
decline, especially after the blistering critique of Lashley and 
Clark (1946). Bailey and von Bonin (1951), for example, 
divided the human neocortex into only 6 areas, even blurring 
the sharp border between VI and V2 originally recognized by 
Meynert (1867). 

Despite subsequent advances in mapping cortical areas in 
non-human primates (Merzenich and Kaas 1980), specifically 
human cortical parcellation languished (though see Braak 
(1980) on "pigmento"-architectonics). The growth of human 
neuroimaging in the 1980s had the effect of elevating Brod- 
mann's map — after a rough translation from Brodmann's 2D 
summary images to a 3D atlas — to a world standard. Despite 
its acknowledged shortcomings (e.g. the middle temporal 
area (MT/V5) did not appear), the moderate resolution of 
early volume-averaged neuroimaging data did not initially 
demand better. Two decades later, as modern postmortem 
human cortical architectonics augmented by immunohisto- 
chemistry began to appear, the situation has improved con- 
siderably (e.g. Malikovic et al. 2007; Caspers et al. 2012). 

However, a persistent stumbling block to combining struc- 
tural and functional data has been to integrate postmortem 
parcellations based on high-resolution microscopic measures 
with retinotopic, connectional, and functional data from in 
vivo brains. Variability between postmortem and in vivo 
brains — whether at the single subject or group level — 
introduces irreducible uncontrolled systematic variation. 

Most previous attempts at in vivo architectonics have only 
been able to consider few localized regions of the cortex 
(Clark et al. 1992; Walters et al. 2003; Eickhoff et al. 2005; 
Sigalovsky et al. 2006; Geyer et al. 2011). One reason is 
signal-to-noise constraints due to in vivo scan time limits. But 
a more serious problem is the difficulty of reliably detecting 
small signal differences (1-4%) between different parts of 
non-quantitative ^-weighted (7\w) images that have large ar- 
tifactual variations in brightness and contrast due to uncor- 
rected inhomogeneities in the local radiofrequency transmit 
and receive fields. Widely applied post hoc histogram-based 
normalization methods (e.g. Dale et al. 1999) and T 1 w/T 2 
weighted (7* 2 w) ratio methods (Glasser and Van Essen 2011) 
remain sensitive to the effects of uncorrected Bl "transmit" 
inhomogeneities (flip-angle variation) on brightness and 
contrast, which affect both within- and between-individual 
comparisons. Ratio methods must also contend with vessel 
artifacts and local spatial distortion that differ between the 2 
scan types, and unfavorable error propagation inherent in 
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Figure 1. Relaxation rate [Hi) as function of cortical depth, area, and curvature, ifl) 
Overall cortical average /?, (equivalent 7", on right /-axis) decreases monotonically 
from gray-white boundary (depth fraction 0.0) to a slight plateau at middle depths 
(0.3-0.6, beginning where second derivative in lower graph crosses zero), and then 
resumes its decrease in superficial layers (0.7-1.0; error bars: ±1 standard error of 
the mean over subjects). (S) Cross-ROI differences in average cortical /?, shown as 
line for 8 depths (from 0.1 near white matter to 0.9 near pia); (/-axes, error bars as 
in ifl); ROIs: Angular (angular gyrus), angular-fs (FreeSurfer angular gyrus label), MT 
low and MT high [non-overlapping low- and high-probability MT labels (Malikovic 
et al. 2007; Fischl et al. 2008), Vl-fs (FreeSurfer V1)]. All matched-paired f-tests on 
hypothesized differences significant (P < 0.05) except where marked "m" [P < 0.1), 
"=" (no significant difference), or "-" (difference opposite prediction). (C) Vertex-wise 
correlation (adjusted R 2 ) of Ri and curvature as function of depth (error bars as 
before over subjects). Scatter plot inset at right is from a single subject at depth 0.5. 
For comparison, left inset shows myelin-stained section of human cortex with reduced 
myelination in concave areas (compare 2 white arrows, from Annese et al. 2004). 

ratio estimation. Finally, differences in local cortical curvature 
(convex vs. concave) systematically affect the laminar and 
myeloarchitecture of the cortex, even within a cortical area 
(Fig. 1; Smart and McSherry 1986; Annese et al. 2004). 

Here, we present a detailed analysis that combines func- 
tional and structural data. To overcome the inherent 



limitations of non-quantitative structural imaging listed above, 
quantitative Tj maps were acquired for all subjects, providing 
a measure of local myeloarchitecture (Schmierer et al. 2004; 
Draganski et al. 2011). The quantification allows for the direct 
comparison of datasets acquired on different subjects at differ- 
ent times and different MRI sites. The Tj maps were then 
combined with retinotopic mapping across multiple sessions 
in the same group of subjects. High-resolution cortical-surface 
reconstructions of the gray-white matter and pial surfaces 
were made from the structural scans and used to analyze, 
visualize, and fuse the 7\ and retinotopy data. Since the longi- 
tudinal relaxation rate, -r?i(=l/7\), positively correlates with 
the level of myelination and the brightness of 7*iW scans as 
commonly used for morphometry, we henceforth refer exclu- 
sively to it. 

Materials and Methods 

Six subjects (ages 22-55, 3 female) with normal or corrected- 
to-normal vision participated in all parts of the study. Experimental 
protocols were approved by local ethics committees, and participants 
gave informed and signed written consent. 

Structural Imaging 

Structural images were acquired on a whole-body Tim Trio system 
(nominal field strength 3T, actual 2.89T, Siemens Healthcare, Erlan- 
gen, Germany), with body transmit and 32-channel receive head coil 
at the Wellcome Trust Centre for Neuroimaging. Proton density- 
weighted (PDw) and TjW images were acquired using an in-house 
multi-echo 3D FLASH pulse sequence (Weiskopf et al. 2011; voxel 
size: 0.8 x 0.8 x 0.81 mm 3 , field of view (FOV) = 256 x 216 x 194 mm, 
matrix = 320 x 270 x 240, repetition time (TR) = 23-7 ms, excitation flip 
angle: 6° (PDw) or 28 (^w)). Acquisition was speeded up by 2x 
GRAPPA parallel imaging in the phase encoding and 6/8 Partial 
Fourier in the partition direction. To improve image quality, 4 gradi- 
ent echoes were acquired (echo delay times (TE) = 2.2, 4.75, 7.3, 9. 85 
ms) after each excitation pulse. Each session consisted of four 10 min 
31 s acquisitions (2 PDw and 2 7\w) and 2 shorter scans to estimate 
field inhomogeneities (see below). Quantitative i?, maps were esti- 
mated from the PDw and 7\w images according to the formalism 
developed by Helms et al. (2008) including a correction for imperfect 
RF spoiling (Preibisch and Deichmann 2009). Recent applications of 
this method have demonstrated its robustness (Helms et al. 2008, 
2009; Draganski et al. 2011). To correct for the effect of radio fre- 
quency (RF) transmit inhomogeneities on i?! maps, maps of the trans- 
mit field Bl + were acquired using a 3D echo-planar imaging (EPI) 
spin-echo (SE)/stimulated echo (STE) method (Lutti et al. 2010, 2012; 
FOV = 256 x 192 x 192 mm 3 , matrix = 64 x 48 x 48, TE SE /TE STE = 39.38/ 
72.62 ms, TR = 500 ms, acquisition time 3 min 48 s), which was cor- 
rected for off-resonance effects using a B0 fieldmap. For further 
details, see Supplementary Information. 

Functional Imaging and Analysis 

Standard r 2 *-weighted EPI scans (3.2-mm isotropic resolution, 128 or 
256 volumes each) and a TjW alignment scan with the same orien- 
tation and slice block center were acquired on a 1.5 whole-body Tim 
Avanto system (Siemens Healthcare) at the Birkbeck/UCL Centre for 
Neuroimaging, with body transmit and 32-channel receive head coil 
(see Supplementary Information for details). Five of 6 subjects com- 
pleted at least 8 retinotopy scans (4 scans in 1 subject) and 4 ipsilat- 
eral mapping scans (2560 volumes per subject) across 3-5 additional 
sessions. See Supplementary Information for details of surface-based 
retinotopy analysis. 

In-house OpenGL/Xlib software drove a video front-projection 
direct-view system that stimulated the visual field to an eccentricity of 
at least 57° of visual angle in all directions from central fixation. 
Balanced (clockwise/counterclockwise, in/out) phase-encoded polar 
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angle (wedge) and eccentricity (ring) stimuli contained 
medium-luminance-contrast colored checkerboards, optical flow 
fields, and simultaneous monitor-for-upside-face-among-right-side-up 
and monitor-for-number-among-letters tasks to maintain continuous 
peripheral attention (eight 64 second cycles). Ipsilateral field 
mapping stimuli used low contrast moving versus stationary gratings 
avoiding the center-of-gaze and the vertical meridian (see Supplemen- 
tary Information for details). 

Cortical-Surface Reconstruction and Sampling ofRj Values 
Within Cortical Ribbon 

Cortical surfaces were reconstructed with FreeSurfer (v5.0.0) from the 
aligned (AFNI 3dAllineate, hand-inspected) average of the 2 high- 
resolution 7\w scans. We initially used quantitative R 1 scans but 
experienced localized segmentation failures because some boundaries 
between the pial surface, the CSF, and the skull have different contrast 
than what is currently assumed by FreeSurfer algorithms (see 
Discussion). 

R 1 data sets were sampled along the normal to each gray-white 
matter surface vertex in steps of 10% of cortical thickness (thickness 
estimated in FreeSurfer; Fischl and Dale 2000) and then smoothed 
tangentially at each depth with a 4-mm full width half maximum 
(FWHM) 2D kernel. The FreeSurfer estimate of local curvature 
(smoothed with a 1-mm FWHM 2D kernel) was used as a linear pre- 
dictor of Ri values at each depth. Vertex-wise residuals from this 
regression were used as "de-curved" estimates of R 1 values whose 
units are directly comparable to raw de-meaned R 1 values. See Sup- 
plementary Information for details of R 1 region of interest (ROI) 
choice. 

Results 

Rj as a Function of Cortical Depth 

Postmortem studies of cortical myelination demonstrated a 
consistent pattern of myelination versus depth across most 
cortical areas. Myelination is highest immediately above the 
white matter, there are 2 bands of high myelination in deeper 
layers (inner/deeper and outer/shallower stria of Baillarger), 
and finally, there is a stepwise reduction in myelination in 
superficial layers. To investigate whether this overall pattern 
was apparent in our R t data, we sampled the R t values for 
each surface vertex at different fractional cortical depths, and 
then averaged these profiles across all vertices. 

Average R± values (Fig. 1A) show an extremely consistent 
decrease from deep to superficial layers (small to large depth 
fractions, z). A moderate plateau in deeper cortical layers 
begins around depth fraction 0.3, where the second spatial 
derivative with respect to z (d 2 R 1 /dz 2 ), shown immediately 
below, crosses zero (an inflection point). The pattern of R t 
with depth closely resembles the overall profile of myelination 
seen in postmortem histology. It is considerably smoother, 
not distinguishing the 2 stria of Baillarger, which is a reflec- 
tion of much coarser sampling units of in vivo human MRI 
(0.8 mm 3 here) compared with histology (~0.01 mm 3 ), as well 
as the averaging of R 1 across the entire cortex of both 
hemispheres. 

Rj Variation over Probabilistically Defined Cortical 
ROIs 

For an initial verification of regional differences in R 1 as a 
measure of myelination, we defined 5 ROIs with large ex- 
pected myelination differences. The first was the FreeSurfer 
probabilistic VI label (Vl-fs). The second and third were MT 
high and MT low (non-overlapping high and low probability 



of cross-subject overlap regions). Finally, we included 2 
angular gyrus labels — a smaller one, angular, defined on our 
subjects (non-visually responsive region just superior to MT) 
and a larger one, angular-fs, from the FreeSurfer parcellation. 

Based on the histological literature, we predicted 1) that R 1 
in VI and the MT ROIs should exceed R ± in the angular gyrus 
ROIs at all cortical depths, and that VI should have the 
highest R 1 , 2) R l in high-probability MT should exceed R 1 in 
surrounding lower probability MT in deeper cortical layers 
(depth fractions 0.1-0.6). The systematic R\ curves across 
areas at different depths illustrated in Figure IB confirmed 
these predictions. All matched-paired ?-tests were significant 
at P<0.05 except where noted, with "m" indicating marginal 
(0.05 <P< 0.10) differences, "=" indicating no significant 
difference, and "-" indicating difference opposite to predicted 
direction. ROI-based analyses showed that at the level of 
macroscopic landmarks, our quantitative estimate of i?i paral- 
lels known regional differences in myelination. 

Relationship ofR t to Local Cortical Curvature 

The human cortex has deep sulci, but also a complex pattern 
of concavity and convexity, including many convex regions 
buried within major sulci. Postmortem studies in humans 
show that myeloarchitecture varies significantly with local 
cortical convexity (Annese et al. 2004); more convex regions 
are thicker and more myelinated, especially in middle and 
upper layers. We therefore examined the relationship 
between R ± and local curvature as a function of cortical 
depth. Figure 1C shows that there is a moderately strong cor- 
relation between local curvature and R 1 at middle cortical 
depth fractions (up to a maximum average adjusted R 2 = 0.14) 
that falls off as one approaches the white matter and the pial 
surface. The bottom right inset scatterplot includes a 
least-squares fit line for R ± against vertex-wise curvature; ver- 
tices in convex regions of the cortex (often but not exclusively 
on gyri) have higher R t than vertices in more concave 
regions. Since i?j varies more at the gray-white matter border 
and at the pial border than it does within the cortex itself, 
small errors in tissue segmentation during surface reconstruc- 
tion can have a large effect on laminar R t estimates near those 
boundaries. The fact that i?i /curvature correlation is strongest 
at middle sampling depths suggests that segmentation errors 
are not responsible for our result. 

Because local curvature is also related to local cortical 
thickness (average vertex-wise adjusted correlation between 
curvature and thickness: Adjusted R 2 = 0.081), local cortical 
thickness might explain some of the variation in R l . However, 
the strong correlation between convexity and R\ survives the 
addition of vertex-wise thickness as a covariate. To minimize 
the effects of local curvature in our subsequent analyses of R\ 
variation across the cortical sheet, we therefore initially re- 
gressed R 1 against curvature, and then used the "de-curved" 
residual values as our estimates of regional differences in cor- 
tical myelination. 

Test-Retest Consistency ofRj Across the Cortex 

The top 3 rows of Figure 2 show medial and lateral views of 
the cortical-surface-based cross-subject average (Fischl et al. 
1999) of i?i from 6 subjects (2 R 1 scans each), displayed on 1 
of the 6 subjects' surfaces. The top row (gray-white matter 
surface), second row (pial surface), and third through sixth 



Cerebral Cortex September 2013, V 23 N 9 2263 




Figure 2. Cross-subject surface average relaxation rate test/retest. Top left 
inset redrawn from Flechsig (1920) shows early myelinating areas in gray. Rows 1-6 
show spherical morph cross-subject average A/?, (difference above mean) projected 
back onto the gray-white, pial, and inflated surfaces of 1 subject, fl, values in each 
subject were sampled at a point halfway between the individual's gray-white matter 
and pial surfaces. Row 3 shows the data before removing variance due to local 
cortical curvature; that variance is removed in all other rows. Rows 5 and 6 show 
excellent scan/rescan reliability of average. Lower middle inset shows same 
magnification lateral occipital cutout with individual subjects' /?, maxima (yellow 
contours), superimposed on MT postmortem probability map in purple (Malikovic 
et al. 2007; Fischl et al. 2008), showing that the individual and average /?, maxima 
differ subtly from the postmortem map. See text for midline and anterior insula 
mask. 



rows (inflated surface) show same magnification views of one 
hemisphere. The heat scale overlay in all rows but row 3 
shows "de-curved" i?i values sampled from voxels lying at a 
cortical depth fraction of 0.5, while row 3 is not "de-curved" 
(but still de-meaned). The effects of "de-curving" can be seen 
by comparing rows 3 and 4 (e.g. the thin line posterior to MT 
disappears). Though regional variations in R t are small 
(<0.04 s _1 ), they are remarkably consistent over repeated 
measurements (bottom 2 rows). R 1 data were masked on the 
midline and in the anterior insula — in the first case because 
single hemisphere surfaces cut into the corpus callosum and 
thalamus, and in second case because of the difficulty of 



accurately reconstructing the gray-white matter surface super- 
ficial to the claustrum. Small regions of high cortical i?j in 
posterior cingulate and anterior insula may thus have been 
omitted here. 

These data show a strong resemblance to Flechsig's (Flech- 
sig 1920) survey of perinatal infant myelogenesis, re-rendered 
in the top left inset of Figure 2. Flechsig found that early mye- 
linating cortical regions (dark shading) tended to be the most 
densely stained for myelin in the adult cortex. These regions 
are clearly visible in our R ± data and include primary motor 
(M-I) and somatosensory (S-I) areas, primary and secondary 
auditory areas, as well as a number of visual areas described 
in more detail in the next section. 



Relation ofRj to Retinotopic Maps 

Visual areas were mapped in the same set of 6 subjects using 
a wide-field direct-view front-projection system. In 2-3 
additional scan sessions for each subject, we mapped polar 
angle, eccentricity, and ipsilateral low-contrast visual field rep- 
resentations (Huk et al. 2002). We then performed 
cortical-surface-based cross-subject averages of complex- 
valued (amplitude and phase) data (Sereno and Huang 2006; 
Hagler et al. 2007) and calculated visual field sign from the 
polar angle and eccentricity averages to identify borders in 
early areas (Sereno et al. 1995). A detailed pattern of corre- 
spondences between and retinotopic borders emerged 
when the 2 data sets were fused on the same surface. 

The retinotopic and R\ averages from both hemispheres in 
Figures 3 and 4 are illustrated in identical lateral, posterior, 
medial, and inferior poses (inflated surface) to aid compari- 
son and are described together. lso-R 1 borders (AR 1 = 0.020 
s _1 ) around R t maxima were traced and superimposed on the 
polar angle data using thin white dashes to aid comparisons 
(see Supplementary Information for a movie that smoothly 
varies view and color scale contrast, more clearly illustrating 
the systematic fine level differences in the data shown in 
Fig. 4). Retinotopic borders (vertical meridian -» circles, hori- 
zontal meridian -* thick black dashes) are illustrated on both 
figures as is the posterior boundary of ipsilateral responses to 
low contrast moving gratings (thick yellow dashes in top row). 

VI (see medial views) was characterized by high R 1 (Affi 
~0.031 s _1 ), and its retinotopically defined borders in both 
hemispheres corresponded almost exactly with a sharp drop 
in R t . Superior to VI, V6 was bilaterally identified by a charac- 
teristic upper-to-lower field transition (moving superiorly and 
posteriorly) on the posterior bank of the parieto-occipital 
sulcus (Pitzalis et al. 2006) with R 1 values similar to VI (Affj 
~0.035 s -1 ) with a shallow R 1 valley between it and VI (see 
contrast ramp movie in Supplementary Information). Interest- 
ingly, Flechsig (1920) identified an early myelinating medial 
area in an almost identical location (Fig. 2, top inset, medial 
view). There was an additional shallower maximum of myeli- 
nation (AR 1 ~0.012 s _1 ) just anterior and superior to V6, con- 
taining mostly lower visual fields that we have tentatively 
labeled V6A (see Gamberini et al. 2011, for finer 
subdivisions). 

In lateral occipital cortex, there was a prominent oval of 
high Ri (Affj ~0.029 s -1 ). We initially expected this might cor- 
respond to MT/V5. However, both polar angle retinotopy (re- 
versal at an upper visual field vertical meridian) and the 
posterior bound of ipsilateral responses were situated near 
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Figure 3. Cross-subject surface average retinotopic maps. Spherical morph average 
polar angle maps are projected back to the inflated right and left hemispheres of one 
subject, and shown in lateral (top), posterior (middle-top), medial (middle-bottom), 
and inferior (bottom) views. The posterior boundary of ipsilateral visual responses is 
marked by a thick yellow dashed line. Thin dotted lines indicate the boundaries of 
regions with high quantitative R : values traced from Figure 4. Vertical and horizontal 
meridians traced from field sign calculations (not shown) are shown as lines of small 
circles and thick black dashes. Generalized visual and multisensory area names (V1, 
V2, V3, VP/V3v, V6, VBA, V8/V01, V3A, MT/V5, FST, LIP (multiple), VIP (multiple), 
PrCu [pre-cuneus visual areaj) were drawn judiciously from the existing conflicting 
literature. 



Figure 4. Cross-subject surface average quantitative H\ (=1/7i) maps. Spherical 
morph average maps of quantitative /?, values sampled at 50% of cortical thickness 
are projected back to the same subject as in Figure 3, and posed and annotated 
identically. Quantitative /?, is illustrated as variation from the mean (A/?,) after 
removing variance due to local cortical curvature. A histogram of A/?, values from 3 
right hemisphere visual areas (VB, MT and VIP defined by retinotopy) and 2 right 
hemisphere non-visual areas (inferior pre-cuneus "default mode" area and non-visual 
lateral occipital cortex) is shown in the lower middle inset. The non-visual regions of 
interest are shown as dashed lines in matching blue and green (upper left and middle 
left). The maxima shown are 3-4% higher than the average /?,. 



the middle of the oval. This region was located within a 
sulcus (Supplementary Fig. 2, showing individual subject data 
on each folded surface; the average was also resampled to 
each subject and shown as a green outline). We labeled the 
posterior part MT. To avoid name proliferation, we labeled 
the anterior and inferior part of this area, fundus of the 
superior temporal sulcus area (FST), on the basis of a similar 
retinotopic layout in non-human primates [e.g. macaque 
monkey, owl monkeys; see also Kolster et al. (2010) and 
Amano et al. (2009) who have subdivided this region into a 
larger number of areas with somewhat different names]. 

In dorsolateral occipital cortex near the midline, we ident- 
ified a patch of high i?i (AR 1 ~0.028 s" 1 ) that corresponded 
well with V3A (defined by the first appearance of upper field 
here moving anteriorly) in the left hemisphere and roughly so 



in the right. Inferiorly, beginning with the center-of-gaze rep- 
resentation of VI (see inferior view), there is a large anteriorly 
extending region of high i?i values (Affi ~0.031 s -1 ) that also 
covers the adjoining center-of-gaze representations V2, 
V3/VP, and hV4 in both hemispheres. This maximum extends 
across parts of at least 4 different visual areas and may be arti- 
factual (see Discussion). Still on the inferior surface, moving 
medially along the horizontal meridian of hV4 (Wandell et al. 
2007), an anterior right-angle bend is visible in a region orig- 
inally identified as V8 by Hadjikhani et al. (1998). This pre- 
cisely corresponds to a small maximum of R\ (Affi ~0.021 
s _1 ) in both hemispheres. 

Back on the lateral surface, there is a prominent maximum 
of i? x (Aff x ~0.020 s _1 ) in both hemispheres at the expected 
location of the ventral intraparietal area (VIP) (Sereno and 
Huang 2006). This is connected to the large high R t strip in 
somatosensory and motor cortex via a small isthmus. In our 
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polar angle mapping, this region contained lower visual field 
and horizontal meridian responses. The lower field bias may 
partly reflect head-centered remapping of receptive fields as a 
result of the lowered gaze position common with direct-view 
mapping. 

Just posterior to VIP, best visible in the posterior view, 
there is an elongated region of moderately high R 1 values 
(AR 1 ~0.013 s _1 ) in both hemispheres that extends posteriorly 
through the region originally identified as the human intra- 
parietal area (LIP) (Sereno et al. 2001) eventually joining up 
with the prominent V3A maximum. Subsequently, this region 
has been subdivided into a number of LIP's (or IPS's) and 
contains several somewhat variable visual maps (e.g. Swisher 
et al. 2007; Silver and Kastner 2009). Finally, in frontal cortex, 
there are shallow R 1 maxima in the frontal eye fields, in an 
area identified as the polysensory zone (Graziano and Gandhi 
2000; Sereno and Huang 2006; Huang and Sereno 2007), 
which appears as an extension off the motor strip, and there 
is a shallow R 1 maximum in the retinotopic part of the dorso- 
lateral prefrontal cortex (Hagler and Sereno 2006). 

Discussion 

Cortical-surface-based analysis of quantitative longitudinal 
relaxation rate, i?i (=1/2"!), showed reliable small (1-4%) 
regional differences in cortical gray matter likely due to differ- 
ences in myelination, as verified by test-retest and by statisti- 
cal analysis in selected regions defined in previous studies 
(postmortem-defined VI and MT, angular gyrus). We com- 
pared this in vivo architectonic measure with retinotopic 
mapping data obtained in separate sessions from the same 
subjects. This revealed that a number of boundaries visible in 
the Ri maps correspond to previously recognized retinotopic 
borders, demonstrating the feasibility of using R 1 maps to 
help with cortical parcellation. Though not all adjoining areas 
have distinguishable R\ values, a number can now be ident- 
ified structurally without elaborate functional scans. 

The cortex-wide average of i?j as a function of laminar 
depth resembled postmortem myeloarchitectonic measures, 
systematically decreasing from the gray-white matter border 
to the pial surface with a moderate plateau in middle 
layers. These results are consistent with findings from loca- 
lized high-resolution imaging and postmortem imaging 
(Walters et al. 2003; Eickhoff et al. 2005). The R t ranking 
of different ROIs — from most to least myelinated — is also 
remarkably consistent over sampling depths. 

R ± was significantly related to local curvature at intermedi- 
ate depths — with convex parts of the cortex more strongly 
myelinated — quantitatively confirming observations in 
myelin-stained sections (Smart and McSherry 1986; Annese 
et al. 2004). Local curvature predicted R^ better than a 
measure of whether a cortical region was on a sulcus or gyrus 
(the primary measure used for cross-subject surface align- 
ment; Fischl et al. 1999). The strength and generality of this 
finding across the cortex suggests that there is a fundamental 
asymmetry between concave and convex bends in the cortex 
having to do with differences in the way deeper and more 
superficial layers can be deformed (Xu et al. 2010). 

At a macroscopic level, the pattern of i?i across the cortex 
showed a remarkable resemblance to postmortem studies of 
the order of myelogenesis (Flechsig 1920). By comparing these 
R 1 maps with retinotopic maps, we found that the borders of a 



number of visual areas defined in the same subjects — including 
VI, V3A, V6, V6A, V8, VIP, retinotopic dorsolateral prefrontal 
cortex, and the frontal eye fields — were marked by abrupt 
increases in R t (0.02-0.04 s _1 ). Perhaps our most surprising 
finding concerned the heavily myelinated oval in the lateral 
occipital cortex. Previously identified as MT/V5 in a number of 
studies, retinotopic mapping (polar angle, posterior boundary 
of ipsilateral responses) showed that MT proper only 
accounted for the posterior 1/3 to 1/2 of that oval. This 
suggests that previous in vivo and postmortem studies in 
humans may have overestimated the extent of MT. Given that 
there are several areas with above average myelination anterior 
and inferior to MT in monkeys (MST and FST — e.g. Bock et al. 
2009), this observation is not unprecedented; however, those 
other areas may be relatively larger or contain additional subdi- 
visions in humans. Somatomotor and auditory cortex were also 
prominently recognizable (treated elsewhere). 

Potential Applications 

In contrast to Tjw imaging, quantitative R 1 imaging directly 
estimates a basic physical property, the longitudinal relaxation 
rate of protons in a given tissue. Thus, the resulting maps are 
directly comparable without intensity normalization across 
individuals, scanners, and time — which especially lends itself 
to studying cortical areal differences, development, health 
versus disease, and disease progression. In this respect, quan- 
titative neuroimaging is superior to standard postmortem 
myelin stains, which, despite their much higher resolution, 
are more subject to the uncontrollable vagaries of silver 
impregnation. 

As an example, R t values in middle cortical layers for the 
high-probability MT label differ significantly from immedi- 
ately adjacent regions, with an average R^ value of ~0.62 s" 1 
at a fractional cortical depth 0.5 (Fig. IB), and an average de- 
curved AR 1 of 0.023 s _1 . To investigate whether this highly 
myelinated region could be identified on a single subject 
basis, we thresholded each subject's surface-smoothed (4-mm 
FWHM) de-curved R 1 map at the average Aff t value of 0.023 
s _1 and searched for a disconnected supra-threshold R\ patch 
within a lateral occipital MT+ search space based on previous 
human studies (Annese et al. 2005). Such patches were 
identified in all 6 subjects, and in 11 of 12 hemispheres 
(Supplementary Fig. 1). 

Methodological Issues and Prospects 

Automated in vivo mapping of architectonics in single sub- 
jects requires detecting local changes in R ± on the order of 1% 
at a high isotropic resolution of ~800 um, which is made 
possible by advanced quantitative methods. In particular, we 
employed high-quality mapping of the B1+ transmit field 
(Lutti et al. 2010, 2012) and highly sensitive multi-echo 3D 
FLASH R t mapping to achieve maximal resolution, accuracy, 
and precision in the shortest time possible (Weiskopf et al. 
2011). We optimized the TR, flip angle and RF spoiling phase 
increment for highest accuracy and precision (Helms et al. 
2011) and applied corrections for imperfect spoiling of mag- 
netization coherence pathways (Preibisch and Deichmann 
2009). No post hoc brightness intensity normalization was 
used. The excellent scan-rescan reproducibility (Fig. 2) 
demonstrates that our R 1 estimates are robust, showing that 
quantitative in vivo maps can be obtained in only 25 min of 
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scan time (see Supplementary Information for further discus- 
sion and description of i?i mapping). 

In recent work done in parallel with ours, cortical myelina- 
tion has been estimated using a ratio between 7*iW and r 2 w 
scans (Glasser and Van Essen 2011). The results were broadly 
similar to ours but focused on large group results. Also, 
an additional series of post hoc normalization steps were 
required to visualize this non-quantitative data, which contain 
uncorrected artifacts in tissue contrast and brightness due to 
non-uniform Bl + transmit fields and to the effects of dividing 
pixel values from the 2 different image types that have con- 
trasts sensitive to different aspects of microanatomy as well as 
different artifacts. The normalization steps are more difficult 
to reproduce at other centers than would be a hard threshold 
at cortical depth = 0.5 of R ± = 0.62 s" 1 . 

Although local increases in cortical R 1 corresponded well 
with known patterns of myelination and with retinotopically 
defined areas, there were some unexpected local increases in 
R 1 — for example, at the nearly adjoined center-of-gaze 
regions of VI, V2, V3, and hV4. Because R 1 differences 
between cortical areas are comparable with R 1 differences 
between different cortical depths (Fig. 1), detection of interar- 
eal differences requires uniformly accurate estimation of the 
gray-white and pial surfaces just as much as it requires i?i 
uniformity. To improve those estimates, the more uniform R l 
images rather than TjW images could be used to reconstruct 
the pial surface. However, this will require extending the 
current FreeSurfer pial-surface-finding algorithm, which is 
finely adapted to 7\w image intensity relations among tissue 
types. For example, in quantitative R ± images, there is no 
change in image intensity between the superior margin of the 
cerebellum and inferior margin of the inferior occipital lobe — 
which leads to a ballooning of the pial surface there — 
whereas in 7\w images, the cerebellum and inferior occipital 
lobe are separated by a thin but distinct hypointensity corre- 
sponding to the dura — likely a mixed contrast artifact, but 
one that FreeSurfer relies on. Another region where the Free- 
Surfer pial-surface reconstruction using R t images fails is near 
the thin temporal bone. It will be necessary to strategically 
incorporate quantitative proton density images into the pipe- 
line to resolve this surface-reconstruction problem. 

Another problem arises in areas such as S-I and VI where 
parts of the cortical ribbon are so thin (~1.5 mm) that they 
are spanned by less than two 0.8 mm 3 voxels. The problem of 
misestimating the cortical depth fraction and hence misesti- 
mating Ri is particularly acute in this case, which might 
explain the unexpected reduction in R t in parafoveal VI (see 
medial surface views in Fig. 3). 

Quantitative i?i mapping is fitted to answering a number of 
questions beyond those introduced here. For example, we are 
currently using it to characterize auditory and somatomotor 
fields (Dick et al. 2011) and are combining it with 
high-angular-resolution diffusion measurements in the cortex. 
This method could also illuminate structure/function 
mapping within cortical areas (e.g. in V2 stripes). Finally, the 
significant individual differences in local quantitative R 1 we 
have uncovered may have functional correlates. 

Supplementary Material 

Supplementary material can be found at: http://www.cercor. 
oxfordjournals.org/. 
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